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Abstract 

The non-compact lattice version of the Abelian Higgs model is 
studied in terms of its topological excitations. The Villain form of 
the partition function is represented as a sum over world-sheets of 
gauge-invariant "vortex" strings. The phase transition of the system 
is then related to the density of these excitations. Through Monte 
Carlo simulations the density of the vortex sheets is shown to be a 
good order parameter for the system. The vortex density essentially 
vanishes in the Higgs phase, and the Coulomb phase consists of a 
single vortex condensate. 



1 Introduction 



In many lattice models the use of topological excitations is helpful in un- 
derstanding the phases of the system. In compact U(l) gauge theory, for 
example, the confining transition is driven by the condensation of magnetic 
monopoles |], [|. These monopoles are consequences of the periodicity in 
the lattice gauge action. Monopoles also appear to be responsible for the 
confinement-deconfinement transition in the compact Abelian Higgs model. 
However, they do not explain the transition separating the Higgs phase from 
the Coulomb and confinement phases |4], [5|, ||. This transition appears to 
be associated with vortex-like excitations, due to the U(l) symmetry in the 
Higgs part of the action [|, |5], [jj . In this paper we will study the role of these 
excitations in the non-compact version of the Abelian Higgs model. Because 
the gauge action is not compact, no monopoles appear in the theory.[] There- 
fore, we are able to show more clearly how the vortex excitations influence 
the Coulomb-Higgs transition. 

The non-compact model is defined as 

/oo poo 
/ D0DA M e- s S = S g + S h (1) 

s b = fE^O*) F^x) = A„A v {x) - A^x) 

X,fll/ 

Sh = -7E(C^^ + M + cc) + EC^ + A5:(C0 x -l) 2 

X,fl X X 

<p x = p(x)e lx ^ x ' is the complex scalar Higgs field, and U X:fl = e lA ^ x \ Sh is the 
same as in the compact model. 

It has become customary in both the compact and non-compact cases to 
consider Sh in the limit where the Higgs self coupling diverges (A — > oo). In 
this limit the magnitude of the Higgs field is constrained to unity, leaving 
only an angular degree of freedom. The resulting total action is 

S = fE^CO - 2 7 Ecos(A mX (x) - A,{x)). (2) 

This fixed-length model is believed to lie in the same universality class, and 
it is simpler to analyze numerically. The phases of the system have already 

1 One can continue to define monopoles exactly as in the compact case. However they 
do not enter into the action and have no bearing on the phases of the system. For a study 
of this see rcf . || . 
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Figure 1: The phase structure of the non-compact Abelian Higgs model, for 
A — > oo. 



been determined in previous Monte Carlo studies [|. A sketch of the phase 
diagram is shown in Fig. [I]. Large values of (3 and 7 correspond to a Higgs 
phase with a massive photon. While small values of (3 or 7 correspond to a 
Coulomb phase, with a massless photon. As f3 — > 00 the model reduces to 
the XY model, which has a continuous transition at j c = 0.15. The 7 — > 00 
limit describes an integer Gaussian model (the "frozen superconductor" of 
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ref. |10|), with a first order transition at f3 c = 0.04. See ref. 
of analytical results. 

Now consider the role of the scalar field in the action. The action must 
be gauge-invariant. Hence, Sh depends upon the phase of the scalar field 
only through the gauge invariant quantity A^x( x ) ~ A^x). Since x(x) is an 
angular variable it can be multi-valued along a "vortex" string (the lattice 
analog of a Nielsen-Olesen string [12[]), and cannot be gauged away. The 
presence of these strings in the theory is very important. Without them 
the fixed-length model would, necessarily, describe just a free gauge boson 
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in the continuum limit. Therefore, these excitations should be central to 
determining the phases of the system at large length scales. 

Einhorn and Savit j|] have shown how to express the partition function 
of the compact model in its dual form, consisting of sums over closed vor- 
tex strings and open vortex strings terminating on monopoles. Polikarpov, 
Wiese, and Zubkov [13] have found a similar expression using the Villain 
form of the non-compact model. In the non-compact case only closed vortex 
strings appear due to the lack of periodicity in S g . In the following section we 
will argue for the existence of a phase transition driven by the condensation 
of these strings, which can be identified with the Coulomb-Higgs transition. 
The intuitive explanation being that the vortices act to disorder the scalar 
field. (As (3 — > oo this is consistent with the results of ref. in four di- 



mensional XY model.) In section 3 the results of numerical simulations will 
be presented demonstrating that this picture is indeed correct. 



2 The Villain Approximation 

The Villain form JTSf of the fixed-length partition function is 

/oo /* oo oo 
/ D%DA M ]T e-^^-TE(AMX-V^) 2 ( 3 ) 

(where we have suppressed the spatial dependence of all the fields). a M is an 
integer- valued field which preserves the periodic symmetry of the original ac- 
tion. This form can be thought of as a "low-temperature" (large 7) Gaussian 
approximation to the action around each minimum of Sh- 

By a series of exact transformations, the partition function above can be 
rewritten as 

00 

Z ~ V" e - 47r2 7l] J^(x)D(x-y;m 2 )J„ l ,(y) ^ 
a M =— 00 

where J^ u = e^xpAxcip is an integer-valued "vortex" current. D(x — y; m 2 ) 
is the usual lattice Green's function, with m 2 = -I. The details of the 
derivation are provided in the appendix. The identification of J^ u as the 
vorticity in the integer part of the Higgs action should be clear. Notice that 
ApJpv = 0. Consequently, J^ u forms closed sheets (world-sheets of closed 
strings) on the dual lattice. 
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Now consider the possible phases of (^) in terms of the vortex excitations. 
Due to the behavior of the lattice Green's function, the self energy of the vor- 
tex sheets will be much larger than any spatially separated interactions. The 
partition function can then be approximated by retaining only the diagonal 
terms. 

oo 

z ~ Y, e- 47r27D(0;m2) £ J ^ (x) (5) 

a M =— oo 

For reasonably large 7, configurations with J lip {x) = 0, +1, —1 will be domi- 
nant. The energy of a vortex sheet of area A will be Ea = 47r 2 7 .D(0 ; m 2 )A. 
Now consider the entropy of such a configuration. Let Na be the number of 
possible closed sheets of area A containing a given plaquette of the lattice. 
For reasonably large A the leading behavior will be Na ~ ^ A -f\ The entropy 
of a sheet of size A will then be Sa ~ Alogfi. Thus, balancing energy and 
entropy produces the critical line 

log// = 47r 2 7 D(0;f), (6) 

where sheets of all sizes are expected to occur .0 

By considering the behavior of D(0 ; m 2 ) it can be seen that the critical 
line above is consistent with the phase diagram of Fig |l]. For example, as 
7 — > 00 D(0 ; m 2 ) ~ ^2 , producing a critical point at 

For larger values of (3 only relatively small surfaces will be present, leading 
to an ordered (Higgs) phase. For smaller (3 values the surfaces are allowed to 
grow large and intersect, forming a single vortex condensate and a disordered 
(Coulomb) phase. The Coulomb-Higgs transition can then be interpreted as 
a condensation of these vortex sheets. 



2 See the second paper of ref. [Q for a derivation of this result, fx > 1. However, a 
good estimate of fi is unknown to the author. 

3 The same critical line exists in the Villain approximation (7, f3 1) to the compact 
model |Q. The two models should agree at large values of (5. In this case, however, no 
approximation has been made to S g so there is no restriction on the values of f3. 
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3 Numerical Results 



Now let's reconsider the original form of the partition function ([!]). The 
conclusions of the previous section can be motivated here by examining the 
periodicity of the lattice Higgs action. For each link of the lattice perform 
the decomposition 

A^x(oc) - A^x) = e^x) + 27ra M (x), (8) 

where 9^{x) G (— tv,tt) and a M (x) is an integer. Notice that changing a^(x) 
on any link of the lattice leaves Sh unchanged. The result is the creation of 
a Dirac string (or vortex) threading through each of the plaquettes of which 
the link is an element. The string current being the same as that defined 
previously^ 

J^{x) = e liV x p A\a p {x). (9) 

If 7 is reasonably large the fluctuations in 9{x) will be small. The entropy of 
the system will be entirely due to J^ u (x). The energy being S g = 2n 2 l3J2 J 2 - 
Thus, we arrive at the same critical energy-entropy balance we found in the 
Villain approximation. 

It is not clear if this can be carried over to smaller values of 7 where the 
fluctuations in 9(x) are not small. To examine the behavior of the vortex 
excitations at these smaller values of 7 Monte Carlo simulations of the fixed- 
length model were performed. The simulations were done on a 8 4 lattice with 
periodic boundary conditions, and in the unitary gauge. The vortex currents 
were then measured directly, using the definition above. 

The simplest measurement which can be made is the amount of vortex 
current per dual plaquette 

v =(w E (10) 

where N p is the number of plaquettes on the lattice. Figures § and |3| show the 
results of such measurements for 7 and (3 ranging from to 1. V essentially 
vanishes in the Higgs phase and rises sharply across the phase boundary to 
a non-zero value in the Coulomb phase. Figure |] shows separate "heating" 
and "cooling" runs at fixed (3 = 0.2. The phase transition, given by the peak 

4 The same definition was used in studying the compact model in ref. |§. Notice the 
similarity with the definition of monopoles in compact QED. In compact gauge theories 
the gauge action is periodic in the plaquette variables. The monopoles are then taken to 
be the oriented sum of the integer parts of the plaquettes around each 3-cube. 
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Figure 2: V plotted as lines of constant vortex density. V = 0.4, 0.3, 0.2, 0.1, 
from left to right. 




Figure 3: Surface plot of the vortex density, V, as measured on a 8 4 lattice. 
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Figure 4: Measurements of V on a 8 4 lattice, at f3 — 0.2. Crosses correspond 
to data taken with increasing 7 steps. Circles correspond to decreasing 7 
steps. 
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in the specific heat, is marked by the arrow at 7 C = 0.24. Statistical errors 
are of the size of the data points. The hysteresis in the data at the critical 
point appears indicative of the continuous (2nd order) phase transition there. 
Runs at other values of the couplings produced similar, excellent, agreement. 
Therefore, V appears to be a good order parameter for the system. 

To better understand the distribution of the vortices, a routine was writ- 
ten to measure the area of each connected vortex sheet. The results support 
the description of the phase transition given in section 2. Deep into the 
Higgs phase only small sheets are present, containing six or ten dual plaque- 
ttes each. As the system is moved towards the transition the average size of 
the clusters increases. Near the transition the clusters begin to coalesce into 
one large cluster, and the number of separate clusters drops dramatically. 
This occurs in the Higgs phase, before the phase boundary is reached. When 
the system reaches the Coulomb phase there is only a single large cluster, 
roughly the size of the lattice. If the system is heated further, the density of 
this cluster continues to increase. 

In conclusion, the numerical simulations presented here seem to confirm 
the picture presented by the large 7 analysis of the theory. In the Higgs 
phase the vortex density is low and the area of each sheet is relatively small. 
In the Coulomb phase the vortex sheets are large and overlap, forming a 
single cluster of dual plaquettes. The abrupt change in the vortex density is 
coincident with the transition, given by the peak in the specific heat. Hence, 
the phase transition in the non-compact Abelian Higgs model, for A — > 00, 
appears to be driven by the condensation of vortex excitations. 

The author would like to thank John B. Kogut for assistance with the 
Monte Carlo simulations and other helpful discussions. The author was sup- 
ported in part by a GAANN fellowship. 



Appendix 

Here it is shown how to transform the Villain form of the partition func- 
tion into equation (f|). The steps are essentially the same as those used in 
analyzing the compact model For greater clarity we will suppress the 
spatial dependence of the fields. (See ref. for a derivation using lattice 
differential forms.) 
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We start with 

Z = D)CD^e"2 E(Wp A A^p) 2 -7E( A MX-^M+27ra M ) 2 ^ ^ 

„ — . — J —oo 



a u =—oo 



and use the identity 



e -2 x =J— dye-^ +iyx (A.2) 



oo 



pa 



to rewrite each of the Gaussian integrals. 

Z ~ f; f^DL^ 

a M =-oo J -°° 

— ^H-i5 pCT e pCTMI ,A MJ 4„-^+iL M (A M x-^M+2™ M )J ^ ^ 

Integrating over x and A p gives 

00 />oo 

z ~ E / DL p DS^A^>(^ + W A ->V) 

The constraint imposed by S(A p L p ) can be removed by letting 

L p = e pupcr A u B pcr . (A. 5) 

The second constraint then becomes 

e^KiBpa + S^) = 0, (A.6) 

implying that 

B pa + S pa = A p ip a - A a ip p . (A.7) 

We still must choose a gauge for B pa to insure that L p is properly defined. 
The most convenient choice is to take B pa + S pa = 0. The partition function 
is then 

00 />oo 

z ~ E / D ^ 

a M =-oo- 7 -°° 
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with m 2 = jj-. Finally, integrating over S pcr , we get the desired result 



Z ~ Z ^2 e" 47 " 27 ^ J l ^(x)D(x-y;rn 2 )J llv 



(y) 



(A.9) 



where (-A 2 + m 2 )D(x - y;m 2 ) = 5 xy and J^ v = e^A^a,,. Z is the 
determinant given by setting all the = in ( |A.8| ), and is completely 
analytic. 
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